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Summary 

A  back-projection  algorithm  is  applied  to  353  Pn  travel  time  measurements  to  image  the 
lateral  variations  in  compressional  velocity  in  the  upper  most  mantle  in  and  around  Tibet  Since 
large-scale  lateral  variations  are  known  to  exist  in  the  Tibetan  crustal  structure,  we  refmed  the 
apparent  Pn  slowness  measurements  prior  to  the  tomographic  inversion.  These  refinements  are 
based  on  a  recently  developed  relocation  procedure  ^ao  &  Helmberger,  1991),  and  on  a  new 
theoretical  understanding  of  head  wave  propagation  in  3D  media  with  velocity  gradients  (Zhao, 
1993).  These  refinements  are  crucial  in  obtaining  a  high-quality  velocity  image. 

We  obtained  an  average  P-velocity  value  for  the  uppermost  mantle  of  7.93  km/s,  and  an 
average  mantle  P-velocity  gradient  of  3.1  X  10*^  srl  for  Tibet  The  final  2D  P-velocity  image  in 
the  uppermost  mantle  is  characterized  by  a  large  (7  to  9  degrees  in  dimension)  low-velocity 
(down  to  about  7.78  km/s)  region  in  north  central  Tibet  and  by  large  high-velocity  (up  to  8.1 1 
km/s)  regions  in  the  western  and  eastern  flanks  of  Tibet 

By  attributing  the  remaining  slowness  residuals  to  systematic  variations  in  crustal  thickness 
among  28  sourc^receiver  regions,  we  estimated  the  crustal  thicknesses  in  those  regions.  The 
result  is  characterized  by  large  thicknesses  (>70  km)  inside  the  Tibetan  plateau,  and  by  decreases 
in  the  crustal  thickness  toward  the  southern,  eastern  and  western  boundaries  of  the  plateau;  the 
thicknesses  at  and  beyond  these  boundaries  are  similar  to  those  in  the  rest  of  Eurasia.  The  rate  of 
the  decreases  appear  to  be  higher  toward  the  western  boundary. 

Rigorous  analyses  of  resolution  and  error  suggest  that 

(i)  there  is  a  pronounced  difference  between  the  tomographic  velocity  image  obtained  with 
and  without  the  refinements  to  the  Pn  slowness  measurements,  indicating  a  significant 
improvement  in  the  quality  of  our  inversion  by  these  refinements; 

(ii)  random  errors  in  measuring  earthquake  origin  times  and  locations,  and  in  the  Pn  arrival 
times,  should  not  cause  errors  of  more  than  0.1  km/s  in  the  velocity  image; 

(iii)  the  resolution  is  generally  within  6  to  10  degrees. 

We  have  also  compared  the  images  in  Pn  velocity  and  in  crustal  thickness  obtained  in  this 
study  with  those  localized  results  from  numerous  pervious  studies  in  and  around  Tibet.  General 
agreement  has  been  found  between  ours  and  most  previous  results.  The  few  exceptional 
disagreements  can  be  interpreted  largely  as  the  result  of  effects  of  a  complex  medium,  such  as 
the  effects  of  mantle  velocity  gradient,  which  have  not  been  carefully  accounted  for  in  previous 
studies. 

The  velocity  image  obtained  in  this  study  is  consistent  with  previous  S-P{  travel  time 
studies  and  geothermal/petrological  studies,  and  is  consistent  with  a  mantle  convection  model  in 
which  a  convective  upwelling  is  occurring  in  north  central  Tibet  and  downwellings  occurring  in 
the  surrounding  regions  (Molnar,  1988, 19^). 
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B^TRODUCTION 

The  Tibetan  Plateau,  the  Himalaya  and  the  Karakoram  are  the  most  spec¬ 
tacular  consequences  of  the  collision  of  the  Indian  subcontinent  with  the  rest  of 
Eurasia  in  Cenozoic  time  (Pig.  1).  Accordingly,  the  deep  structure  beneath  them 
proiddes  constraints  on  both  the  tectonic  hbtory  of  the  region  and  on  the 
dynamic  processes  that  created  these  structures  (Molnar,  1088). 

There  have  been  numerous  studies  of  seismic  velocity  structure  beneath 
Tibet.  Studies  of  surface  waves  (  eg.  Bourjot  &  Romanowicz,  1002), 
waveforms  (  Holt  &  Wallace,  1000),  and  S-P  differential  travel  times  (Molnar  & 
Chen,  1084;  Pettersen  &  Doombos,  1087;  Molnar,  1000)  have  indicated  that  there 
are  rignificant  lateral  variations  in  seismic  velodties  in  and  around  Tibet.  There 
seemed  to  be  a  tendency  for  velocities  to  be  lower  in  north  central  Tibet  (Molnar, 
1000).  This  tendency,  while  inferred  from  very  sparse  data  coverage,  has  been 
intriguing  because  of  its  suggestion  of  a  probable  convective  upwelling  in  north 
eentral  Tibet  and  downwellings  in  the  surrounding  regions  (Molnar,  1088,  1000), 
which  are  in  agreement  with  the  observations  of  Cenozoic  volcanism  in  north  cen¬ 
tral  Tibet  (Burke  et  al.,  1074;  Deng,  1078;  Kidd,  1075;  Sengor  &  IGdd,  1070; 


Molnar  et  a!.,  1087). 

Previous  tomographic  inverrions  of  travel  times  have  been  conducted  in 
the  western  U.S.  (eg.,  Hearn  et  al.,  1001).  The  P  wave  velocity  images  from  these 
inver«(ms  generally  correlate  well  with  the  major  tectonic  features  in  those  areas. 
In  P,  tomography,  it  b  desirable  to  invert  for  the  lateral  variations  of  P  wave 
velodty  in  the  uppermost  mantle.  When  there  b  a  velocity  gradient  in  the  upper 
mantle,  however,  P,  observed  at  different  distances  samples  mantle  P-wave  veloci> 
ties  through  different  depth-ranges.  Thb  has  been  shown  to  be  the  ease  in  Tibet 
(Song  et  al.,  1085;  Holt  &  Wallace,  1000).  Previous  P.  tomography  studies  have 
ignored  the  effect  of  the  mantle  velocity  gradient.  In  an  area  where  thb  gradient 
exists,  images  of  lateral  variations  in  P  wave  velocity  may  be  contaminated  by 
vertical  variations  in  the  mantle  P  wave  velocity.  Recently  Zhao  (1003)  developed 
mathematical  formulations  which  enable  us  to  calculate  head  wave  travel  times  in 
the  presence  of  lateral  structural  variations,  and  to  estimate  and  correct  for  the 
upper  mantle  velocity  gradient,  if  this  gradient  b  nearly  constant  with  the  vary¬ 
ing  depth. 

In  this  study,  we  apply  the  mathematical  formulations  of  Zhao  (1003),  as 
well  as  the  relocation  procedure  by  Zhao  &  Helmberger  (1001)  to  44  seismic  events 
and  353  P,  travel  times  to  obtain  P  wave  slownesses  along  paths  confined  to  the 
uppermost  mantle.  These  slownesses  are  then  input  to  a  tomographic  back- 
projection  algorithm  to  obtain  a  laterally  varying,  2D  P  wave  velocity  image  in 
the  uppermost  mantle  beneath  the  Tibetan  Plateau.  Under  a  further  assumption 
that  the  remnant  slowness  residual,  which  b  not  diminished  by  tomography,  b 
caused  primarily  by  the  vertical  travel  time  delays  at  the  sources  and  stations,  we 
also  invert  for  the  crustal  thicknesses  for  28  sub-regions  in  Tibet.  We  also  con¬ 
duct  rigorous  analysis  of  error  and  resolution  of  the  P  wave  velocity  inversion, 
and  test  the  sensitivity  of  the  2D  P  velocity  image  to  (i)  event  relocation,  (ii)  the 


CMTcctloDs  of  effects  of  large-scale  lateral  structural  variatioss,  and  (Hi)  the  correc¬ 
tions  of  effects  of  mantle  velocity  gradient.  We  will  demonstrate  that  (i)  plays  a 
role  in  correctly  resolving  the  detailed  lateral  variations  in  the  2D  P-velodty 
image,  whereas  (ii)  and  (iii)  mainly  affect  the  averse  (DC  level)  of  the  P  velocity 
image.  Fmally,  the  2D  P  velocity  image  and  the  crustal  thicknesses  will  be  com¬ 
pared  to  the  numerous  localized  studies  conducted  prior  to  this  study,  and  will  be 
discussed  in  terms  of  tectonic  processes  underneath  the  Tibetan  Plateau. 

DATA 

We  collected  353  P,  arrival  times  from  44  earthquakes  (Table  1)  which 
occurred  inside  Tibet.  Fig.  2  shows  the  ray  coverage.  The  distance  range  of  the 
data  is  3-15  degrees.  P„  arrivals  at  greater  distances  were  not  used,  since,  from 
synthetics  using  model  TIP,  which  was  developed  for  Tibet  by  Zhao  et  ai.  (1901), 
deep  turning  P  waves  arrive  before  P„  beyond  18  degrees.  Of  these  travel  times, 
we  repicked  42  to  judge  the  reliability  of  the  ISC  (International  Seismological 
Centre)  picks.  Generally,  they  agree  within  0.5  seconds.  We  rejected  the  P„ 
arrival  times  yielding  apparent  velocities  less  than  7  km/s  or  greater  than  0 
km/s. 

The  parameters  of  the  events  used  in  this  study  are  obtained  by  the  reloca¬ 
tion  procedure  of  Zhao  &  Helmberger  (1001).  These  parameters  are  listed  in 
Table  1.  The  number  of  the  P«  phase  picks  from  the  event  is  also  listed  in  Table 
1,  under  In  Table  2,  we  listed  the  stations  used  in  thb  study,  with  vertical 
travel  time  delays  {V  j)-  The  latter  are  calculated  from  the  available  regionalized 
models,  including  the  models  of  Holt  &  Wallace  (1000),  Zhao  et  al.  (1001),  and 
Zhao  &  Helmberger  (1003).  More  details  of  the  relevant  calculation  will  be  dis¬ 
cussed  later. 


INVERSE  METHODS 

Computerized  tomography  has  previously  been  applied  to  travel  times  to 
image  the  P  velodty  In  the  upper  mantle  (eg.,  Hearn  et  al.,  1001).  In  general,  tfj, 
the  P^  travel  time  from  an  Ith  source  to  a  ith  station  can  be  expressed  by 

*ij  -  +  oi+tj  ,  (!) 

where  Xfj  b  the  source  to  receiver  horizontal  distance,  Vy  b  the  apparent  P  wave 
velodty,  a,-  and  6y  are  vertical  travel  time  delays  at  the  source  and  receiver, 
respectively  ^eam  et  al.,  1001).  Zhao  (1003)  has  shown  that  this  relation  is  valid 
for  the  head  waves  traveling  in  complex  waveguides,  except  that  a,-  and  bj  should 
be  calculated  using  local  structures. 

Method  to  estimate  a,-,  6y  and  mantle  velocity  gradient 

In  previous  studies,  a,-  and  bj  are  solved  for  iteratively.  For  this  study,  a 
friori  knowledge  of  the  averaged  velocity  structures  within  and  outside  Tibet  is 
available  (Holt  &  Wallace,  1000;  Zhao  et  al.,  1091;  Zhao  &  Helmberger,  1003). 
Our  approach  b  therefore  to  use  this  a  priori  knowledge  to  estimate  a,*  and  bj 
approximately,  and  subtract  these  approximate  estimates  from  the  travel  times 
before  the  tomographic  inversion,  s,*  and  bj  estimates  thus  obtained  do  not 
include  the  effects  of  smaller^scale  lateral  variations  in  crustal  thickness,  since  we 
do  not  have  a  priori  knowledge  of  these  variations. 

Another  difference  between  our  method  and  previous  inverse  methods  is  that 
in  the  current  inversion,  we  take  into  account  the  P  velocity  gradient  in  the  upper 
mantle,  which  b  known  to  exbt  beneath  the  Ubetan  Plateau  (Song  et  al.,  1085; 
Holt  &  Wallace,  1000).  An  important  effect  of  the  mantle  velocity  gradient  is 
that  the  P.  phase  observed  at  different  distances  samples  P  velocities  in  different 
depth-ranges,  and  the  resulting  P  velocity  image  from  tomography  will  be 


somewhat  contaminated  by  the  depth  variations  of  P  velocity  in  the  mantle. 


Numerically,  both  the  approximate  corrections  for  a,-  and  6j-  and  the  correc¬ 
tions  for  mantle  gradient  are  possible  only  with  the  new  formulations  recently 
devdoped  by  Zhao  (1003).  These  formulations  related  the  P,  travel  times  to  the 
F  vdodties  in  the  uppermost  mantle  in  the  presence  of  both  lateral  structural 
variations  and  a  constant  velodty  gradient  in  the  upper  mantle.  To  use  these  for¬ 
mulations,  we  rewrite  equation  (1)  as 


1  ,  <0  -^1 
V-  *  X- 

•'ij 


(2) 


Denoting  the  estimates  of  a,-  and  bj  by  o',-  and  V  respectively,  we  have 

Vij  ^  0  ^0 


(3) 


When  a  constant  velocity  gradient  is  present,  the  mantle  P-velocity  distribution 
can  be  expressed  as 


V(r)- V(zo)(l+e(r -*o))  »  (4) 

where  V(z)  b  the  depth-varying  P  velocity,  Zq  b  the  depth  of  Moho,  and  c  b  a 
constant  whose  product  wirh  V(ro)  gives  the  mantle  velocity  gradient  for  a  flat 
Earth.  For  a  spherical  Earth,  the  velocity  gradient  b  (c  —  1.6X10“^)  V(ro)  when 
the  units  for  e  and  V(zo)  are  km“^  and  km/z^  respectively  (c/.,  Helmberger, 
1073).  Zhao  (1003)  shows  that  when  such  a  constant  vdocity  gradient  exists, 
in  equations  (1)  and  (2)  can  be  expressed  as 

where  b  the  P  vdocity  in  the  uppermost  mantle  averaged  along  the  {ij)  th 
path,  and  X  ij  b  the  horizontal  distance  that  a  ray  travtU  in  the  mantle.  Now  we 
introduce  a  small  residual 


m  -  (vS-v(io))(i+c»jr^/«) 

» (•) 

'where  V’(«o)  is  VfJ  aversged  over  sll  of  the  («',/)  psirs.  Substituting  equation  (6) 
into  (5)  leads  to 

Va  -  v(i,)  P) 

In  practice  we  must  use  V  (equation  (3))  instead  of  V^j  when  equation  (7)  is 
used;  will  in  that  ease  contain,  in  addition  to  the  right-hand  ride  of  equation 
(6),  contributions  from  non-zero  Sa^  and  6b jt  defined  as 


Sci  -  s'  ,  -e, 

(8) 

1 

1 

• 

(») 

For  an  area  where  a  number  of  different  V  fj  are  collected,  V{sq)  and  e  can  be 
solved  by  linear  regression  using  equation  (7). 

Method  used  to  obtain  2D  P  velodties  in  the  uppermost  mantle 

Once  F(so)  and  c  are  solved,  can  be  estimated  for  each  of  the  (t,y)  pairs 
using  equations  (3)  and  (5).  Dividing  the  area  under  study  into  cells,  we  have  from 
the  definition  of  Vjj 

+  6ie9ij  ,  (10) 

ViJ  ^  a  t 

where  is  the  line  segment  over  which  the  (t,y)th  path  overlaps  the  kih  cell, 
§it  is  the  slowness  in  the  kih  cell,  ReSij  is  the  slowness  residual,  which  is  caused 
by  o'rors  in  event  origin  times/locations,  in  the  travel  time  readings,  and  in  s',- 
and  b^  j,  which  are  approximations  of  a,-  and  bj  and  are  used  to  calculate  VfJ^ 


(equation  (3)).  Strictly  apeaking,  deviations  of  the  real  uppermost  mantle  form  a 
gridded,  smoothly  varying  surface  will  also  contribute  to  e,y;  effects  of  this  kind 
are  known  as  tb^^  "modellsatlon  error"  (eg.,  Tarantola,  1087).  Equation  (10)  forms 
a  standard  finear  inverse  problem  for  unknown  s^.  We  use  a  variation  of  the 
back-projection,  or  SIRT  algorithm  (eg.,  Humphry  Sl  Clayton,  1088),  described 
*ln  detail  by  Xle  &  Mitchell  (1000)  to  solve  for  In  particular,  we  use  the  0- 
polnt  smoothing  operation  (Suetsugu  &  Nakanbhl,  1085;  Xle  Si  Mitchell,  1000) 
between  the  Iterations.  This  operation  Is  used  to  overcome  Insufficient  spatial  data 
coverage,  at  a  slight  cost  of  resolution. 


Method  to  estimate  crustal  thicknesses 

In  the  back-projection  calculation,  values  are  updated  during  the  iterations 

to  reduce  the  misfit  between  the  Inverse  of  and  — — SAi»***  When  the 

^  ij  k 

misfit  stops  decreasing,  we  end  up  with  the  final  slowness  model  (final  values  of 

$i)  and  some  non-diminished  Retij  values.  To  the  first  order  approximation  we 

* 

assume  that  Res^j  values  are  dominately  caused  by  o',-  and  V  j,  the  impredse 
estimates  on  a,-  and  bj  prior  to  the  back-projection  inversion,  we  then  have 


Re$ij 


Sa^  +  Sbj 


(11) 


At  the  time  we  obtain  Reo^j,  a,-  and  bj  are  known,  and  o,-  and  (therefore  5a,- 
and  Sbj)  are  not.  Nbn-sero  5e,-  and  Sbj  values  primarily  result  from  the  lack  of 
smaller-scale  variations  In  the  crustal  structure  used  to  calculate  ef and  V  j  prior 
to  the  inversion.  With  the  known  ReOfj  values,  we  can  solve  for  5a,-  and  Sbj  to 
obtain  Of  and  bj  values,  which  can  then  be  used  to  infer  the  smaller  scale  crustal 
thicknesses  (this  is  just  a  calculation  inverse  to  one  with  which  we  obtained  o'  ,- 
and  V  j).  To  make  the  inversion  for  5a,-  and  Sbj  overdetermined,  and  to  reduce 


ibe  probable  effects  of  errors  in  origin  times  (Table  1),  in  event  locations  (Zhao  & 
Hdmberger,  1001)  and  in  P,  arrival  time  readings  (which  also  contribute  to 
merge  the  unknowxis  £a,-  and  Sbj  into  a  new  set  of  unknowns.  To  do 
so  we  divide  the  area  under  study  into  sub-regions  with  the  following  properties: 
(1)  each  of  the  /th  sub-region  contains  either  a  group  of  sources  or  a  group  of  sta- 
,  tions,  called  a  "sub-region  of  sources'*  or  a  "sub-region  of  receivers";  (2)  within 
each  sub-region,  all  of  the  ef  ,•  and  j  calculated  for  locations  inside  it  have  been 
constant,  and  (3)  the  geographic  extent  of  each  of  the  /th  sub-region  is  small,  so 
that  the  crustal  thickness  can  be  assumed  constant.  We  now  assume  that  Sa,  and 
S6j-  are  constants  within  a  given  sub-region,  le., 

Sa,  -  Ti  (12) 

for  all  of  the  fth  sources  located  within  the  /th  sub-region  of  sources,  or 

Sb^  -  T,  (13) 

for  all  of  the  yth  station  located  within  the  /th  sub-region  of  stations,  f/*  thus 
introduced  form  a  new  set  of  variables  which  is  smaller  in  size  than  the  unknowns 
&,■  and  Sbj,  provided  that  the  number  of  sub-regions  is  small.  Using  RtSij,  we 
can  now  solve  for  r/  in  an  overdetermined  least-squares  problem. 

AVERAGE  Pn  VELOCITY  AND  MANTLE  VELOCITY  GRADIENT 

The  d  i  values  are  obtained  with  a  laterally  homogeneous  velocity  model 
(model  TIP  of  Zhao  et  al.,  1001).  Values  of  y  are  calculated  from  the  models  of 
Holt  &  Wallace  (1000)  and  of  Zhao  et  al.  (1001).  J  ,•  and  V  j  are  all  estimated 
using  the  method  of  Zhao  (1003).  The  resulting  V  j  values  are  listed  in  Table  2. 

The  of  i  and  V  j  values  thus  obtained  are  then  subtracted  from  tij  values  to 
obtain  apparent  P„  velocity  values,  which  were  subsequently  multiplied  by  a  fac¬ 
tor  of  {Re~P)1^e  ^  correct  for  a  first-order  effect  of  the  earth’s  sphericity. 


where  (>^6371  km)  and  D  (wSS  km)  are  defined  as  the  average  radius  of  the 
earth  and  the  average  crustal  thicknesses  in  the  area  under  study,  respectively. 
The  resulting  apparent  P,  velocity  in  equation  (3))  values  are  plotted  in 
Fig.  3  against  distance  {X  fj).  In  Fig.  3,  the  trend  of  increasing  V  fj  with  X  ,-y  is 
obvious,  indicating  the  existence  of  the  mantle  P  velocity  gradient.  We  used  the 
method  described  in  the  last  section  to  invert  for  V(so)  and  c,  (Eq.  (7)).  The 
results  are 


V(io)- 7.03  ±  0.17km /s  , 

(H) 

e  -  5.5  X  10“^ifcm“* 

(IS) 

For  X  ij  less  than  about  1700  km,  it  is  clear  that  the  condition  in  equation  (6)  is 
valid.  The  V(ro)  and  e  values  thus  estimated  predict  the  solid  line  in  Fig.  3, 
with  the  dashed  lines  representing  the  standard  deviations.  Equation  (14)  gives, 
for  the  uppermost  mantle  in  Tibet,  an  averaged  P  velocity  that  agrees  well  with 
7.04  km/s  by  Song  et  al.  (1085).  This  value  b,  on  the  other  hand,  lower  than 
8.11  km/s  by  Jia  et  al.  (1081)  and  8.12  km/s  by  Chen  &  Molnar  (1081).  The 
diflference  may  be  caused  by  the  relocation  and  correction  of  mantle  velocity  gra¬ 
dient,  both  are  conducted  in  this  study.  Our  value  of  7.03  km/s  is  abo  lower 
than  8.24  km/s  by  Holt  &  Wallace  (1000)  and  8.20  km/s  by  Zhao  et  al.  (1001), 
and  both  were  obtained  by  using  long-period  body  waveforms  which  sample 
deeper  parts  of  the  mantle. 

Equations  (14)  imd  (15)  mean  that  the  average  velocity  gradient  in  the  litho¬ 
sphere  beneath  the  Tibetan  Plateau  and  its  surrounding  regions  b  3.1X10~^s~* 
tor  the  upper  150  km  of  mantle,  after  the  effect  of  the  Earth’s  sphericity  b 
corrected.  Song  et  al.  (1085)  obtained  a  depth-dependent  velocity  model  for  the 
upper  mantle  in  Tibet.  We  have  applied  a  least  squares  fitting  to  their  velocity 


valxioi  in  the  upper  150  km  of  mantle  and  obtained  a  velocity  gradient  of 
2.0X10~*«*’*,  which  agrees  well  with  our  estimate.  Holt  &  Wallace  (1090) 
obtained  a  rough  estimate  of  2.5X10“*s“*  for  the  same  velocity  gradient,  a  value 
not  significantly  diflerent  from  our  estimate. 

•  LATERAL  VARIATIONS  IN  P,  VELOCITIES  AND  CRUSTAL 
THICKNESS 

The  c  value  in  equation  (15)  and  sf if  j  terms  obtained  were  used  to  esti¬ 
mate  values  (equation  (5)),  which  were  then  used  in  the  back-projection 
tomography  to  solve  for  (equation  (10)).  We  divided  the  area  in  and  around 
Tibet  covered  by  the  ray  paths  in  Fig.  2  into  984  cells  of  a  size  of  1.5*  by  1.5*. 
Convergence  in  the  back-projection  calculation  takes  place  after  20  iterations. 
The  resulting  2D  P  velocity  image  is  shown  in  Fig.  4. 

To  estimate  the  lateral  variations  in  the  cnistal  thickness  using  the  inverse 
method  described  earlier,  we  divided  the  area  under  study  into  28  sub-regions. 
We  set  a  linear  least-squares  problem  to  solve  for  r/,  which  are  then  used  to  esti¬ 
mate  crustal  thicknesses  within  the  sub-regions  (Table  3  and  Fig.  5). 

The  overall  pattern  of  the  2D  P  velocity  image  in  Fig.  4  is  characterized  by  a 
low  velodty  (down  to  about  7.78-7.83  km/s)  region  in  north  central  Tibet,  and 
high  velocity  regions  in  the  north-eastern  and  western  flanks  of  Tibet.  This  corre¬ 
lates  well  with  the  regional  variations  in  S-P  travel  times  reported  by  Molnar  & 
Chen  (1984),  Pettersen  &  Doombos  (1987)  and  Molnar  (1990).  In  Fig.  5,  many  of 
the  areas  inside  Tibet  have  great  (>  70  km)  crustal  thicknesses.  Crustal  thick¬ 
ness  generally  decreases  toward  and  beyond  the  boundaries  of  the  plateau;  the 
rate  of  decrease  appears  to  be  lower  toward  the  eastern  boundary  than  toward  the 
other  boundaries.  Next,  we  compare  regional  values  of  our  P  velocity  image  and 
crustal  thicknesses  with  the  previous  localized  studies. 


For  western  Tibet  (west  to  80*E),  including  Hindu  Kush  and  Pamir,  we 
obtained  P  wave  velodties  in  the  uppermost  mantle  of  between  about  8.03  and 
8.11  km/see  figure  4),  and  the  crustal  thicknesses  of  between  about  65  and  70 
km  (Figure  5).  These  results  agree  well  with  previous  values  obtained  within  the 
same  area  using  various  methods  (Beloussov  et  al.,  1080;  Kaila,  1081;  Holt  &  Wal> 
lace,  1000;  Zhao  St  Helmberger,  1003).  Beneath  station  KSH  we  obtained  a  eru> 
stal  thickness  of  44  km  (Group  11  of  Table  8;  Figure  5).  This  station  is  located 
near  the  western  margin  of  the  Tarim  Basin,  where  the  crustal  thickness  changes 
rapidly  QPeng,  1085;  Zhao  et  al.,  1001).  This  probably  explains  why  Feng  et  al. 
(1081)  and  Feng  (1085)  suggested  a  greater  crustal  thickness  of  50  km  inside  the 
Basin. 

Along  the  foothills  of  the  Himalaya  ( the  southern  boundary  of  the  Tibetan 
Plateau),  we  obtained  P,  velocities  of  7.85  •  8.01  km/sec.  These  values  agree  well 
with  the  results  from  travel  time  studies  by  Kaila  at  al.  (1068),  Verma  (1074)  and 
that  from  waveform  inversion  by  Holt  St  Wallace  (1000).  Our  value  is  also 
consistent  with  the  values  obtained  for  the  Indian  subcontinent  from  deep  seismic 
sounding  (DSS)  studies  (eg.,  Kaila,  1086;  Kaila  et  al.,  1000a,  1000b).  Beyond  the 
southern  boundary  of  Tibet,  we  obtained  cnistal  thicknesses  of  34  ±  4  km  and 
38  ±  4km  for  two  western  locations,  and  35  ±  3  km  and  20^:3  km  for  two  central 
locations  (Fig.  5).  The  latter  values  agree  well  with  the  estimated  thickness  by 
Kaila  et  al.  (1068)  and  Verma  (1074).  Our  estimates  of  crustal  thickness  also 
agree  well  with  those  by  Holt  St  Wallace  (1000)  with  an  exception  near  the 
western  end  of  the  southern  boundary,  where  the  thickness  by  Holt  St  Wallace 
(1000)  is  about  8  km  greater  than  ours. 

In  centra]  Tibet  (80*-00*E),  we  obtained  low  P,  velocities  (7.70  to  7.05 
km/s).  Holt  St  Wallace  (1000,  Fig.  Ob)  suggested  a  much  higher  velocity  (8.24 
km/s).  Their  high  velocity  may  result  from  the  effect  of  velocity  gradient  in  the 


mantle.  Our  estimated  erustal  thicknesses  ranse  between  72  and  80  km,  eon* 
aistent  with  71.0  km  by  Holt  Sc  Wallace  (1000). 

In  the  Lhasa  block,  there  have  been  three  DSS  profiles,  l.e.,  1  (28.8*N, 
8S.5*E)  -  (28.8  "N,  00.3^),  2  (32.8*N,  01.4*E)  •  (27.4*N,  80.2^),  and  3  (32.0*N, 
80.1*E)  •  (32.(fN,  04.5*E),  during  the  Sino*Franoo  Joint  Project.  For  this  block, 
we  obtained  vdodties  of  7.88  to  8.00  km/s,  and  crustal  thicknesses  of  75  km 
in  the  north,  and  50  km  in  the  south.  The  erustal  thickness  from  interpreting  the 
DSS  profiles  were  70*77  km  along  profile  1  (eg..  Him  et  al.,  1084a;  Teng  et  al., 
1083b,  1085b;  Xiong  et  al.,  1085);  75  to  40  km  from  north  to  south  along  profile  2 
(eg..  Him  et  al.,  1084b,  1084c;  Him  Sc  Sapin,  1084;  Institute  of  Geophysics,  1081; 
Teng  et  al.,  1883a);  and  70  km  along  profile  3  (Sapin  et  al.,  1085;  Teng  et  al., 
1085b).  These  agree  well  with  our  estimates.  Along  profile  1,  Teng  et  al.  (1083b, 
1085a)  suggested  P,  velocities  of  8.1*8.2  km/sec,  whereas  Him  et  al.  (1084a)  sug* 
gested  a  much  higher  value  of  8.7  km/s.  Along  profile  2,  Teng  et  a\  (1081; 
1083a}  and  Institute  of  Geophysics  (1081)  estimated  P„  velocities  of  7.05*8.4 
km/s.  These  wide  ranges  in  the  estimates  of  P,  velocities  maybe,  due  to  many 
factors,  including  the  ambiguity  of  P,  arrivals,  the  dipping  Moho,  error  in  the 
erustal  velocity  stractures  used,  and  the  narrow  distance  range  involved. 

In  southeastern  Tibet  (south  to  35*N  and  east  to  05*£),  including  the  Yun¬ 
nan  and  Sichuan  provinces,  we  obtained  P„  velocities  that  increase  from  7.81 
km/s  in  the  southeast  to  7.05  km/s  in  the  northwest.  The  crustal  thickness 
estimated  for  the  interior  of  this  area  is  subject  to  fairly  large  uncertainties  (40  ± 
17  km;  see  Fig.  5).  This  b  because  only  7  arrivals  are  available  for  the  estimate 
0;roup  21  in  Table  3).  The  real  averaged  erustal  thinkness  near  that  location  may 
be  well  above  60  km.  Our  estimated  crustal  thickness  is  40  ±  5  km  beneath  sta¬ 
tion  KMI;  and  is  51  ±  6  km  beneath  stations  CDN  and  CD2  (Fig.  5).  These 
results  agree  with,  within  the  uncertainties,  the  P,  velocities  and  crustal 


tbleknesBcs  by  Holt  &  Wallace  (1000).  Previous  studies  have  suggested  a  changing 
crustal  thickness  from  37  to  70  km,  as  well  as  changing  velodty  from  7.75  to 
8.10  km/sec  (eg.,  Feng,  1085;  Hu  et  al.,  1086;  Zhu  et  al.,  1086,  Xiong  et  al.,  1086; 
Zhao  &  Zhang,  1087;  Yin  et  al.,  1088;  Kan  &  Un,  1088)  for  this  area. 

To  summarize  the  above  comparisons,  the  distributions  of  both  F,  velocity 
and  crustal  thickness  in  Figs  4  and  5  generally  ^ree  with,  within  uncertainties, 
previous  results  wherever  available.  It  is  interesting  that  the  low  ve  s  in 
north  central  Tibet  coincide  not  only  with  the  S>P  residual  lows  previously 
obtained,  but  also  with  the  locations  where  young,  Cenozoic  volcanic  rocks  have 
been  observed  (Burke  et  al.,  1074;  Kidd,  1075;  Deng,  1078;  Sengor  &  Kidd,  1070; 
Molnar  et  al.,  1087).  These  support  the  hypothesis  by  Molnar  (1088,  1000)  that  ^ 
convective  upwelling  may  be  occurring  in  north  central  Tibet,  and  downwelling 
may  be  occurring  in  surrotmding  regions. 

ANALYSIS  OF  RESOLUTION  AND  ERROR 

In  this  section  we  quantitatively  assess  the  quality  of  the  2D  P^yelocity  image 
obtained  in  this  study.  We  also  analyze,  in  a  less  rigorous  manner,  the  uncer^ 
tainty  of  the  crustal  thicknesses  reported  in  the  earlier  section. 

Resolution  of  the  P  velocity  image 

In  computerized  tomography,  the  spatial  resolution  at  a  given  cell  can  be 
approximated  by  the  point  spreading  function  (p.s.f.)  calculated  for  the  cell  (eg., 
Humphreys  &  Clayton,  1088;  Xie  &  Mitchell,  1000).  We  calculated  the  p.s.f  for 
many  cells  in  Tibet.  Figs  6a,  6b  and  6c  show  p.s.f.  at  three  representative  cells 
located  in  western,  central,  and  eastern  Tibet,  respectively.  As  illustrated  by 
these  examples,  the  spreading  of  the  p.s.f  (and  therefore  the  resolution)  is  gen¬ 
erally  within  about  10*.  At  some  locations  where  the  ray  coverage  is  relatively 


denser,  the  pjbS  npproaehes  to  leas  than  6*  (Fig.  6b). 


Effects  of  random  errors  on  P  vdodiy  image 

The  random  errors  in  event  ioeations  and  origin  times,  and  in  the  P.  travd 
time  readings  can  cause  random  errors  in  the  ^  P  velodty  image.  The  events 
used  in  this  study  have  been  relocated,  the  errors  in  origin  times  are  given  in 
Table  1.  Since  there  is  a  direct  relationship  between  the  errors  in  event  origin 
times  and  the  errors  in  evoat  locations,  the  errors  listed  in  Table  1  actually 
represent  combined  errors  In  all  of  the  event  parameters,  including  both  event  ori¬ 
gin  times  and  in  event  locations  (Zhao  &  Helmberger,  1001).  There  are  also  prob¬ 
able  errors  in  the  ISC  picks  of  the  P,  arrival  times,  for  which  we  estimated  a 
maximum  amplitude  of  1  s.  The  estimated  maximum  amplitude  of  the  probable 
errors  in  our  own  picks  is  0.5  s.  We  think  these  are  good  estimates  since  the 
events  are  fdriy  large  (m)  >  5.5).  To  estimate  the  random  errors  in  the  P- 
velodty  image  caused  by  these  errors  and  uncertainties,  we  used  the  estimated 
amplitudes  for  the  latter  and  a  pseudo>random  sign  generator  (Xie  &  Mitchell, 
1000)  to  generate  synthetic  errors  in  travel  time  measurements.  For  details  and 
advantages  of  using  this  random  error  generator,  the  reader  is  referred  to  Xle  & 
Mitchell  (1000). 

Figure  7  shows  the  resulting  estimates  of  the  random  errors.  For  areas  inside 
Tibet,  the  estimated  random  errors  are  less  than  about  0.02  km/s.  The  errors 
increase  toward  the  boundaries  of  Tibet  because  of  poorer  ray  coverages,  but  the 
values  of  errors  there  are  still  quite  small  (less  than  0.05  km/s),  with  two  excep¬ 
tional  larger  values  in  the  northwest  and  southwest  comer.  We  note,  however, 
that  the  estimated  random  errors  might  be  smaller  than  the  "modelisation  errors'* 
including  those  mentioned  before,  and  those  caused  by  the  probable  anisotropy. 


Sendtivity  of  2D  P  vdodty  image  to  mantle  velodty  gradient 


To  aee  bow  aensltive  the  2D  P-velodty  image  in  Fig.  4  b  to  the  correction  of 
mantle  vdodty  gradient,  we  conducted  back-projection  inversion  using  apparent 
Fs  vdodties  (V*  in  equation  (3))  without  correcting  for  that  gradient,  i.  e.,  we 
replaced  V^j  in  equation  (10)  by  V  fj  in  the  sensativity  test.  The  resulting  P- 
vdodty  image  b  plotted  in  Fig.  8.  In  tbb  figure,  the  main  features  of  the  lateral 
variations  in  P-vdodty  image  b  roughly  the  same  as  that  from  inversion  with  the 
correction  for  the  gradient  (Fig.  4).  However  the  average  velocity  in  Fig.  8  b 
higher  (8.04  ±  0.18  km/s).  Thb  vdodty  b  dose  to  those  obtained  by  Jia  et  al. 
(1081)  and  by  Chen  &  Molnar  (1081),  but  lower  than  those  by  Holt  &  Wallace 
(1000)  and  that  of  8.20  km/s  by  Zhao  et  al.  (1001)  who  used  long-period  waves. 
From  Figures  4  and  8,  corrections  of  the  mantle  vdodty  gradient  do  not  greatly 
aflect  the  patterns  of  lateral  variations  In  the  tomographic  images  of  F„  velocity, 
but  affect  significantly  the  average  values  of  the  images. 

The  upper  mantle  vdodty  gradient  may  have  lateral  variations,  which  can 
not  be  resolved  from  the  current  data  set.  Based  on  P  wave  travd  time  tomogra¬ 
phy  for  volcanic  regions,  Zhao  et  al.  (1002)  suggested  that  an  upwelling  may  have 
a  higher  velocity  gradient  by  about  1.0X10~’s~*  than  the  surrounding  regions.  If 
thb  b  true  for  north  central  Tibet,  the  vdodties  beneath  the  region  would  be 
even  lower  than  in  Fig.  4. 

Effects  of  vertical  travel  times  at  sources  and  recdvers 

As  described  earlier,  we  used  model  TIP  and  the  laterally  varying  velocity 
modeb  by  Holt  &  Wallace  (1000)  to  obtain  ,•  and  V  y,  the  estimates  of  vertical 
travd  times  from  sources  and  stations  to  Moho.  These  estimates  were  then  used 
to  estimate  apparent  P,  velocities,  V  ;y,  which  in  turn  are  used  in  the  back- 
projection  inversion.  To  see  the  effect  of  these  and  V  j  values  on  the  final  P 


▼dodty  imase,  we  eondueted  a  test  oi  inversion  in  which  we  estimated  tif  ,•  and 
It  j  using  model  TIP  only,  which  results  in  a  constant  value  of  6.9  s  for  vertical 
travel  times  from  surface  to  Moho.  The  mantle  velocity  gradient  te  not  corrected 
in  this  testing  invention,  which  results  in  the  P  velodty  image  plotted  in  Fig.  9. 

The  average  velocity  di  8.19  ±  0.17  km/sec  in  Fig.  9  b  even  higher  than  that 
in  Fig.  8.  The  main  feature  of  a  slower  velodty  region  in  northern  central  Tibet 
in  Figs  4  and  8  seems  to  be  preserved  in  Fig.  9,  but  the  velocities  near  the  boun- 
daiy  of  Tibet  in  Fig.  9  becomes  faster.  Thb  b  particularly  pronounced  for  celb 
near  the  western  end  and  the  middle  of  the  southern  boundary.  This  phenomenon 
can  be  explained  by  a  thiner  crust  near  the  boundaries  of  Tibet,  where  the  crustal 
thicknesses  can  be  as  small  as  30  to  35  km.  Model  TIP  predicted  a  universal  70 
km-thick  crust,  which  is  most  improper  near  the  boundaries  and  lead  to  over* 
estimation  of  mantle  velocities  there.  Nevertheless,  the  effects  of  ignoring  the 
mantle  velocity  gradient  and  the  lateral  variations  in  vertical  8tation>to>Moho 
travel  times  have  not  significantly  altered  the  main  pattern  of  lateral  variations  in 
P  velocity  image,  and  the  average  value  of  the  P  velocity  image  seems  to  be 
affected  much  more  profoundly. 

Effect  of  event  relocation  on  P  velocity  image 

In  thb  study,  event  origin  times  and  locations  are  obtmned  by  the  relocation 
procedure  of  Zhao  &  Helmberger  (1991).  Previous  P,  studies  used  event  parame¬ 
ters  and  P,  arrival  times  directly  form  the  ISC  bulletins.  Thb  has  been  the  ease 
for  both  P,  tomography  (eg.,  Hearn  et  al.,  1991)  and  for  the  studies  in  and 
around  Tibet  (Chen  &  Molnar,  1981;  Barazangi  &  Ni,  1982;  Ni  8i  Barazangi, 
1983).  Zhao  and  Helmberger  (1991)  found  that  the  locations  by  ISC  lead  to  an 
over-estimation  of  about  1.5%  for  velocity  in  Tibet.  Holt  &  Wallace  (1990) 
abo  found  that  ISC  locations  of  Tibetan  events  were  systematically  too  deep,  thus 


kmdlng  to  liigher  velodty  estimates.  To  test  the  effects  of  not  conducting  event 
rdocation,  we  conducted  the  back-projection  inversion  using  the  parameters  by 
ISC  bulletins,  with  corrections  of  V j  calculated  using  the  model  TIP,  and 
without  the  corrections  of  mantle  velocity  gradient.  The  resulting  P  velocity 
image  is  plotted  in  Fig.  10,  where  the  average  velocity  is  8.28  ±  0.10  km/s, 
which  is  nmilar  to  8.42  km/s  suggested  by  Barazangi  St  Ni  (1082)  and  by  Ni  & 
Barazangi  (1083).  The  main  pattern  of  a  slower  north  central  Tibet  is  no  longer 
dear  in  Fig.  10,  and  many  of  the  faster  velodtles  in  the  flanks  of  Tibet  in  Figs  4, 
8  and  0  are  now  outside  Tibet.  Fig.  10  gives  one  a  misconception  that  nearly  all 
of  Tibet  forms  a  large  low-velocity  region.  This  test  therefore  indicates  that  it  is 
extremely  important  to  relocate  the  events  in  a  P,  tomography  study,  at  least  for 
an  area  like  Tibet. 

Errors  in  the  crustal  thickness  estimates 

The  estimates  of  the  crustal  thickness  plotted  in  Fig.  5  are  obtained  under 
the  assumption  that  the  remnant  travel-time  residual,  after  the  backrprojection,  is 
dominantly  caused  by  Soi  and  Sbj  in  equations  (8)  and  (0),  which  reflect  the 
small-scale  variations  in  crustal  thickness  near  the  source  and  stations.  In  making 
such  an  assumption,  we  have  ignored  the  probable  effects  on  the  remnant  residu¬ 
als,  lystematicaliy  caused  in  each  of  the  28  sub-region,  by  (i)  origin  time  errors, 
Oi)  erroneous  P,  arrival  readings,  and  (iii)  the  deviation  of  the  real  Moho  inter¬ 
face  from  the  gridded,  smoothly  varying  surface.  These  assumptions  cannot  be 
perfectly  valid.  However,  the  fact  that  high-level  structural  variations  have  been 
known  to  exist  in  Tibet  (eg..  Him  &  Sapin,  1084)  suggests  that  our  assumptions 
diould  be,  to  the  flrst  order,  valid  and  that  the  crustal  thickness  thus  obtained 
should  serve  as  first-order  approximations. 


Another  practical  Imperfection  In  the  result  of  of  our  Inversion  for  the  crustal 
thicknesses  Is  that  the  result  may  not  be  reliable  for  groups  with  a  small  number 
of  picks  (under  '!N*'  in  Table  3).  For  one  point  In  Fig.  S  (near  40*N,  0S*E),  the 
crustal  thickness  Is  estimated  to  be  >23  db  17  km,  resulting  from  only  7  P, 
travel>tlme  redduals.  We  have  found  that  the  corresponding  Inverse  matrix  Is 
ringular.  We  think  that  the  crustal  thicknesses  Inferred  from  other  groups  with 
less  than  10  arrivals  (Table  3)  may  not  be  reliable. 

DISCUSSION  AND  CONCLUSIONS 

353  P„  travel  times  from  44  Tibetan  earthquakes  (Table  1)  and  at  46  seismic 
stations  (Table  2)  have  been  collected  to  study  lateral  variations  in  compresslonal 
velocity  in  the  uppermost  mantle  underneath  the  Tibetan  plateau.  The  inverse 
method  and  procedures  of  this  study  differ  from  those  in  previous  P.  tomography 
studies  in  that  in  the  current  study  we  applied,  to  the  travel  time  measurements, 
corrections  to  several  important  biases,  such  as  those  caused  by  (i)  event  mislocs' 
tion  by  ISC,  Cii)  mantle  velocity  gradient  and  (Hi)  large-scale  variations  in  crustal 
thicknesses. 

We  used  a  back-projection  algorithm  to  obtain  a  2D  P  velocity  image  for  the 
uppermost  mantle.  We  also  inverted  for,  to  the  first  order  approximation,  the 
crustal  thicknesses  in  28  sub-regions  in  and  around  Tibet. 

The  main  results  of  the  Inversion  can  be  summarized  in  the  following  para¬ 
graphs: 

(1)  The  average  P  velocity  value  for  the  uppermost  mantle  In  Tibet  is 
7.03  ±  0.17  km/s. 

(2)  The  average  P-velocity  gradient  in  the  upper  150  km  of  the  mantle  Is  3.1  X 
10“®  8“*. 


(3)  The  2D  P'velodty  image  of  the  region  includes  a  low-velodty  tone  in  north 
central  Tibet,  and  two  high*velodty  sones  in  the  western  and  eastern  flanks 
of  Tibet. 

(4)  In  much  of  the  area  inside  Ubet,  the  crustal  thickness  is  greater  than  70  km. 
Near  and  outside  the  boundaries  of  the  plateau,  there  is  a  general  trend  of 
decreasing  crustal  thickness,  and  the  rate  of  decrease  appears  to  be  lower 
toward  the  eastern  boundary  than  toward  the  other  boundaries. 

Comparison  of  our  2D  P-velocity  image  and  estimates  of  crustal  thicknesses 
with  those  obtained  in  previous  localized  studies  shows  general  agreement  within 
the  uncertainties.  A  few  exceptions  can  be  expiained  by  various  complications  in 
propagation,  such  as  the  effects  of  mantle  velocity  gradient,  which  has  not 
been  taken  in  to  account  in  most  previous  studies. 

We  conducted  rigorous  analyses  of  the  resolution  and  random  error  in  the  2D 
P  velodty  image.  The  spatial  resolution  of  that  image  is  estimated  to  be  gen¬ 
erally  within  10®,  and  often  within  6®  where  the  ray  coverage  is  dense.  The  ran¬ 
dom  error  in  that  image  is  estimated  to  be  genially  within  about  0.02  km/s  for 
locations  inside  Tibet.  However,  the  systematic,  'modelisation"  errors  might  be 
greater.  Our  estimates  of  the  crustal  thicknesses  are  reliable  only  to  the  first-order 
because  of  the  assumptions  used,  and  the  often  limited  number  of  slowness  residu¬ 
als  available  in  a  sub-region. 

Another  important  finding  of  this  study  is  that  event  relocation  (correction 
to  bias  0)  mentioned  above)  plays  a  very  important  role  in  reliably  retrieving  the 
detailed  lateral  variations  In  2D  P  velocity  image.  Corrections  to  biases  (ii)  and 
Oil)  mentioned  above,  on  the  other  hand,  seem  to  have  a  greater  effect  on  the 
average  (DC  level)  of  the  velocity  image;  the  average  P  velocity  is  over-estimated 
without  these  corrections. 


The  low-vdodty  region  in  north  eentn)  Tibet  found  in  this  study  coincides 
with  the  locations  of  low  S-P  travel  time  residuals  found  in  localised  studies  by 
Molnar  and  Chen  (1084),  Pettersen  and  Doombos  (1087)  and  Molnar  (1000).  It 
abo  overlaps  with  the  area  where  Neogene  or  younger  volcanic  rocks  have  been 
found  ^urke  et  al.,  1074;  Kidd,  1075;  Deng,  1078;  Sengor  &  Kidd,  1070;  Molnar 
ct  a!.,  1087).  These  support  the  hypothesis  by  Molnar  (1088,  1000)  that  a  convec* 
tive  upwelling  may  be  occurring  in  north  central  Tibet,  and  downwelHng  occur* 
ring  in  surrounding  regions. 
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Table  1.  Earthquakes  used  in  this  studyf 


Date  Origin  Time  Location  Depth 


Event 

AfoPa  Yf 

hr  min 

aee 

(•N 

•F) 

(km) 

Nt 

1 

3 

16  64 

1 

5 

13.5  ±  1.2 

37.11 

95.60 

10 

2 

2 

3 

6  66 

2 

10 

52.1  ±  0.8 

31.50 

80.53 

11 

4 

3 

3 

6  66 

2 

15 

49.8  ±  1.1 

31.48 

80.50 

8 

3 

4 

10 

14  66 

1 

4 

40.0  ±  1.2 

36.50 

87.46 

8 

7 

5 

8 

15  67 

9 

21 

•2.6  ±  0.9 

31.23 

93.56 

8 

2 

6 

8 

30  67 

4 

22 

01.8  ±  1.2 

31.71 100.24 

10 

1 

7 

8 

30  67 

11 

9 

44.8  ±  1.4 

31.73  100.28 

8 

1 

8 

3 

24  71 

13 

54 

15.2  ±  1.2 

35.46 

98.03 

7 

1 

9 

4 

3  71 

4 

49 

•2.2  ±  1.4 

32.19 

95.08 

8 

1 

10 

5 

22  71 

20 

3 

27.3  ±  1.3 

32.42 

92.11 

8 

3 

11 

8 

30  72 

15 

14 

05.0  ±1.3 

36.64 

96.35 

15 

1 

12 

8 

30  72 

18 

47 

38.0  ±  1.2 

36.71 

96.46 

19 

2 

13 

2 

6 

73 

10 

37 

06.1  ±  1.5 

31.38100.54 

10 

1 

14 

2 

7 

73 

16 

6 

20.7  ±  1.3 

31.58100.28 

8 

1 

15 

7 

14 

73 

4 

51 

16.1  ±  1.5 

35.25 

86.41 

6 

12 

16 

7 

14 

73 

13 

39 

24.1  ±  0.7 

35.30 

86.51 

7 

10 

17 

11 

18 

77 

5 

20 

06.7  ±  1.1 

32.75 

88.44 

10 

8 

18 

3 

29 

79 

7 

7 

15.3  ±  1.3 

32.55 

97.28 

10 

7 

19 

2 

22 

80 

3 

2 

42.8  ±  0.7 

30.65 

88.61 

10 

9 

20 

6 

1 

80 

6 

19 

•5.2  ±  1.1 

39.06 

95.63 

12 

7 

21 

6 

24 

80 

7 

35 

44.7  ±  1.2 

33.02 

88.46 

11 

10 

22 

7 

29 

80 

14 

58 

39.4  ±  1.0 

29.64 

81.08 

15 

7 

23 

8 

23 

80 

21 

36 

48.9  ±  1.5 

33.07 

75.69 

14 

6 

24 

8 

23 

80 

21 

50 

•0.5  ±  2.0 

32.89 

75.83 

13 

6 

25 

10 

7 

80 

9 

32 

03.0  ±  1.3 

35.63 

82.21 

4 

7 

26 

11 

19 

80 

19 

0 

45.5  ±  0.9 

27.47 

88.80 

14 

8 

27 

1 

23 

81 

21 

13 

46.5  ±  1.2 

30.99  101.14 

8 

5 

28 

6 

9 

81 

22 

8 

17.6  ±  1.3 

34.65 

91.34 

10 

9 

29 

1 

23 

82 

17 

37 

25.8  ±  1.2 

31.74 

82.27 

10 

4 

30 

6 

15 

82 

23 

24 

26.5  ±  1.2 

31.91 

99.92 

7 

11 

31 

2 

13 

83 

1 

40 

08.2  ±  1.9 

40.05 

75.17 

5 

11 

32 

5 

20 

85 

15 

11 

35.6  ±  0.5 

35.62 

87.28 

8 

20 

33 

6 

20 

86 

17 

12 

43.5  ±  1.2 

31.23 

86.79 

15 

13 

34 

7 

6 

86 

19 

24 

20.9  ±  1.2 

34.55 

80.15 

5 

22 

35 

8 

20 

86 

21 

23 

51.2  ±  1.2 

34.65 

91.49 

11 

15 

36 

8 

26 

86 

9 

43 

•1.2  ±  1.4 

37.83  101.55 

7 

7 

37 

2 

25 

87 

19 

56 

33.3  ±  0.8 

38.12 

91.13 

15 

15 

38 

12 

22 

87 

0 

16 

38.5  ±  1.0 

41.40 

89.66 

20 

13 

39 

11 

5 

88 

2 

14 

29.7  ±  1.1 

34.34 

91.79 

8 

21 

40 

4 

15 

89 

20 

34 

06.9  ±  0.9 

30.04 

99.13 

8 

15 

41 

4 

25 

89 

2 

13 

19.5  ±  0.9 

30.05 

99.43 

8 

14 

42 

5 

3 

89 

5 

53 

•0.7  ±  0.9 

30.07 

99.49 

9 

12 

43 

5 

3 

89 

15 

41 

30.0  ±  0.7 

30.13 

99.55 

8 

12 

-i4- 

_a. 

_2. 

_2i. 

Ei'jbnR'i 

_12_ 

-J. 

t  The  event  parameters  are  obtained  by  the  relocation  procedure  of  Zhao  &  Helmberger 
(1091).  In  particular,  the  latitudes  and  longitudes  are  determined  using  travel  times,  and 
the  depths  are  determined  by  waveform  modeling  of  p,  pP  and  aP  phases.  For  deteails, 
see  Zhao  &  Helmberger  (1991). 


I  N  is  the  number  of  P,  arrivals  available  from  the  event. 


T%bl«  2.  Stations  used  In  tins  study 


Name 

Location 

TD 

(•N 

•E) 

(mc) 

Nt 

AJM 

S6.467 

74.650 

4A 

8 

BRCl 

r.667 

81.580 

8J 

2 

BMI 

87J00 

92.200 

8A 

8 

BHK 

SI  .417 

76.417 

4A 

4 

CDS 

S0.910 

103.758 

5.1 

10 

CDU 

S0.680 

104.011 

5.1 

6 

CEP 

38.824 

71.609 

4J 

2 

CHA 

28.838 

87.167 

8A 

4 

CHOP 

83.658 

73.260 

4J 

8 

CPA 

82.680 

71.424 

4J 

2 

DDI 

80.817 

78.051 

44 

5 

DMN 

27.600 

85.106 

84 

10 

DRP 

81.744 

70.203 

44 

1 

GAR 

89.000 

70.817 

6.5 

11 

GTA 

86.400 

96.800 

6.0 

11 

GtVH 

26.100 

01.467 

8.6 

5 

GYA 

26.450 

106.664 

8.5 

9 

INR 

27.050 

93.283 

4.4 

7 

jm 

26.733 

04.167 

5.1 

4 

JMU 

82.717 

74.900 

8.6 

2 

KHM 

25.651 

04.068 

5.1 

8 

KHO 

37.483 

71.533 

6.5 

13 

KKN 

27.784 

85.288 

8.6 

10 

KKR 

29.951 

76.817 

44 

6 

KMI 

25.123 

102.740 

4.1 

16 

KOI 

26.983 

95.500 

5.1 

5 

KSH 

39.450 

75.968 

6.0 

15 

LAH 

31.550 

74.333 

44 

8 

LSA 

29.700 

91.150 

6.9 

24 

LZH 

86.084 

103.834 

6.0 

19 

MNL 

33.135 

78.750 

44 

3 

NDl 

28.683 

77.217 

44 

24 

NIL 

33.650 

78.250 

44 

6 

PKI 

27.571 

85.409 

8.6 

9 

PSH 

38.937 

71.434 

44 

1 

PTH 

29.550 

80.217 

54 

4 

SARP 

81.922 

72.670 

4.8 

8 

SBDP 

82.300 

70.800 

4.8 

1 

SHL 

25.567 

91.888 

84 

86 

SRNI 

83.950 

74.750 

44 

1 

THW 

82.794 

71.748 

44 

8 

TOC 

26.750 

04.767 

5.1 

5 

VAR 

25.300 

83.017 

84 

14 

WRS 

84.150 

71.401 

44 

5 

XAN 

84.034 

108.917 

8.6 

18 

YYI 

26.567 

94.683 

5.1 

2 

1 /  is  vertical  travel  time  from  Mobo  to  ctatioo. 
t  N  is  Dumber  of  P,  travd  times  used  to  calculate  t’y 


Table  8.  Vertical  Travel 


Time  Corrections  for  Sub^regions  {r,) 


Sub>  region 

Stations  or  Events 

1 

GAR4<HO 

0.78±0.52 

24 

2 

CPA, CEP, CHOP  J)RP,JMU0.84^.41 

LAHa^,NILJ»SH,SBDP 

SRNI,SARP,THW,WRS 

41 

3 

AJM3HKJ)DI4a<R,NDI 

PTH 

-0.41d:0.37 

46 

4 

BRCI,CHA3KN4<MN3KI-0.70±0.35 

VAR 

40 

5 

BMI,GWH,INT,JHI3HM 

KOI,SHL,TOC,YYI 

-0.13^.30 

70 

6 

GYA3MI 

>0.06±0.50 

25 

7 

CD2,CDU 

O.O2±O.50 

16 

8 

XAN 

-0.36±0.66 

13 

0 

LZH 

0.80±0.54 

10 

10 

GTA 

-0.28±0.71 

11 

11 

KSH 

>1.42±0.61 

15 

12 

ev25,ev34 

0.77±0.45 

20 

13 

ev2,ev3,ev22,ev20 

0.88±0.57 

18 

14 

ev26 

•O.60±O.82 

8 

15 

ev4,cvl5,evl6,ev32 

0.47±0.37 

40 

16 

evl7,evl0,ev21  ,ev33 

O.40±O.38 

40 

17 

ev31 

0.18±0.72 

li 

18 

ev37 

0.06±0.60 

15 

10 

ev20,ev35,ev30 

0.08±0.36 

45 

20 

ev5,ev0,evlO 

O.85±O.05 

6 

21 

evl8 

•1.17±O.80 

7 

22 

ev6,ev7,evl3,evl4,ev27 
ev30,ev40,ev41  ,ev42 
ev43,ev44 

0.08±0.30 

80 

23 

evl,ev8,evll,evl2 

O.28±O.08 

6 

24 

ev20 

-4.40±O.88 

7 

25 

LSA 

0.65±0.48 

24 

26 

ev23,ev24 

-0.68±0.70 

12 

27 

ev36 

0.65±0.88 

7 

28 

ev38 

0.30±0.65 

13 

t  N  U  number  of  travel  time  renduals  (S§'f  or  Si*/)  need  to  estimate  r|. 


Figure  1.  3D  view  of  the  ta^iogniphy  of  the  Tibetan  Plateau  and  its  surrounding  regions. 
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Figure  2.  Ray  paths.  Triangles  are  stations.  Stars  are  earthquakes.  A  heavy  solid  line  is 
the  contour  of 3000  m  altitude,  which  nuirks  the  boundary  of  die  Tibetan  Plateau. 
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Figure  4.  2D  P  vclocUy  image  for  the  uppermost  mantle  underneath  the  Tibetan  Plateau. 
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Map  of  the  estimated  random  errors  in  the  2D  P  velocity  image.  Inside  Tibet^  the 
less  than  0.02  km/sec.  the  error  on  the  Tibetan  boundary  is  generally  less  than 


Figure  8.  2D  P  velocity  image  obtained  without  correcting  for  the  mantle  velocity  gradient 


Figure  9.  2D  P  velocity  image  obtained  with  no  corrections  to  the  mantle  velocity  gradient, 
nor  to  the  regionally  varying  vertical  travel  times  from  source/station  to  Moho. 


